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The phenomenom of emerging regular spectral features from random interactions 
is addressed in the context of the interacting boson model. A mean-field analysis 
links different regions of the parameter space with definite geometric shapes. The 
results provide a clear and transparent interpretation of the high degree of order 
that has been observed before in numerical studies. 



1 Introduction 

Recent shell model calculations for even-even nuclei in the sd shell and the 
pf shell showed, despite the random nature of the two-body matrix elements, 
a remarkable statistical preference for ground states with angular momentum 
L = 1 . A similar dominance of L = ground states was found in an analysis 
of the Interacting Boson Model (IBM) with random interactions 2 . In addition, 
in the IBM there is strong evidence for both vibrational and rotational band 
structures. According to the conventional ideas in the field, the occurrence of 
regular spectral features is due to a very specific form of the interactions. The 
studies with random interactions show that the class of Hamiltonians that lead 
to these ordered patterns is much larger than is usually thought. 

The basic ingredients of the numerical simulations, both for the nuclear 
shell model and for the IBM, are the structure of the model space, the ensemble 
of random Hamiltonians, the order of the interactions (one- and two-body), 
and the global symmetries, i.e. time-reversal, hermiticity and rotation and 
reflection symmetry. The latter three symmetries of the Hamiltonian cannot 
be modified, since we are studying many-body systems whose eigenstates have 
real energies and good angular momentum and parity. It was found that the 
observed spectral order is a rather robust property which does not depend on 
the specific choice of the (two-body) ensemble of random interactions 1,3 ' 4 ' 5 , the 
time-reversal symmetry 3 , or the restriction of the Hamiltonian to one- and two- 
body interactions 6 . This suggests that that an explanation of the origin of the 
observed regular features has to be sought in the many-body dynamics of the 
model space and/or the general statistical properties of random interactions. 

The purpose of this contribution is to investigate the distribution of ground 
state angular momenta for the IBM in a Hartree-Bose mean- field analysis 7 . 
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2 Phase transitions 



The IBM describes low-lying collective excitations in nuclei in terms of a sys- 
tem of N interacting quadrupole (eft) and monopole (s^) bosons 8 . The IBM 
Hamiltonian spans a wide range of collective features which includes vibra- 
tional, rotational and 7 unstable nuclei. The connection with potential energy 
surfaces, geometric shapes and phase transitions can be studied by means of 
Hartree-Bose mean-field methods 9 ' 10 in which the trial wave function is writ- 
ten as a coherent state. For one- and two-body interactions the coherent state 
can be expressed in terms of an axially symmetric condensate 



N,a) = -7= (cosas f +sina<4) |0> , (1) 



/M 

with — 7r/2 < a < n/2. The angle a is related to the deformation parame- 
ters in the intrinsic frame, /3 and j 9 . First we investigate the properties of 
some schematic Hamiltonians that have been used in the study of shape phase 
transitions. 

2.1 The U(5)-SO(6) case 

The transition from vibrational to 7 unstable nuclei can be described by the 
Hamiltonian 

H = gg* dt . d+ m S ™ X _ 1} (sV _ rft . J) (gg _ d-d) , (2) 

which exhibits a second order phase transition at \c = 7r/4 9 . For the present 
application, we extend the range of the angle \ to that of a full period — 7r/2 < 
X < 37r/2, so that all possible combinations of attractive and repulsive inter- 
actions are included. The potential energy surface is given by the expectation 
value of H in the coherent state 

E(a) — cos x sin 2 a + i sin x cos 2 2a . (3) 

The equilibrium configurations are characterized by the value of a = o>o for 
which the energy surface has its minimum. They can be divided into three 
different classes or shape phases 

a = -7r/2 < x < tt/4 

cos2a = cotx tt/4 < x < Stt/A (4) 

a = tt/2 3tt/4 < x < 3?r/2 
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which correspond to an s-boson or spherical condensate, a deformed conden- 
sate, and a <i-boson condensate, respectively. The phase transitions at the 
critical points %c = ?r/4 and 37r/4 are of second order, whereas the one at 37r/2 
is of first order. 

The angular momentum of the ground state can be obtained from the 
rotational structure of the equilibrium configuration, in combination with the 
Thoulcss-Valatin formula for the corresponding moments of inertia 10 . 

• For ao = the equilibrium configuration has spherical symmetry, and hence 
can only have L = 0. 

• For < a < tt/2 the condensate is deformed. The ordering of the rotational 
energy levels L = 0, 2, . . . , 2N 

E m t = ^rHL + 1) , (5) 

is determined by the sign of the moment of inertia 

3A^(sinx-cosx) 

-£-3 = : • (.0) 

sin x cos X 

For 7r/4 < x < the moment of inertia J3 is positive and hence the ground 
state has angular momentum L = 0, whereas for for tt/2 < \ < 37r/4 it is 
negative corresponding to a ground state with L — 2N. 

• For a = tt/2 wc find a condensate of iV quadrupole or d-bosons, which 
corresponds to a quadrupole oscillator with N quanta. Its rotational structure 
is characterized by the labels r, n\ and L. The boson seniority r is given by 
t = 3riA + A = N, N — 2, . . . , 1 or for N odd or even, and the values of the 
angular momenta are L = A, A + 1, . . . , 2A — 2, 2A 8 . In general, the rotational 
excitation energies depend on two moments of inertia 

E Iot = ±-t(t + 3) + ±L(L + 1) . (7) 

For the special case of the Hamiltonian of Eq. (2) only the first term is present 

2N 



smx 



(8) 



For 37r/4 < x < tt the moment of inertia X5 is negative and the ground state 
has t — N, whereas for tt < x < 37r/2 it is positive and the ground state has 
t = (L = 0) for N even, and r = 1 (L = 2) for N odd. 

In Fig. 1 we compare the percentages of ground states with L = and 
L = 2 as a function of N obtained exactly (solid lines) and in the mean-field 
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analysis (dashed lines). The results were obtained by assuming a constant 
probability distribution for x on the interval — n/2 < \ < 3tt/2. We have 
added a small attractive L ■ L interaction to remove the degeneracy of the 
ground state for the r — N solution. There is a perfect agreement for all 
values of N. The ground state is most likely to have angular momentum 
L = 0: in 75% of the cases for N even and in 50% for N odd. In 25% of 
the cases, the ground state has the maximum value of the angular momentum 
L = 2N. The only other value that occurs is L = 2 in 25% of the cases for 
N odd. The oscillation in the L — and L = 2 percentages is due to the 
contribution of the d-boson condensate. The sum of the L = and L = 2 
percentages is constant (75%) and does not depend on N. 

2.2 The U(5)-SU(3) case 

A second transitional region of interest is the one between vibrational and 
rotational nuclei. In the IBM, it can be described schematically by 



smx 
V - 1 

+ (2s t xd t ±V7d t xd t )^-(2JxS±V7Jx d)^ . (9) 



H± = cos X J ■ d+ [(2s t -s t -d t -d t )(2S-S- d- d) 
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Figure 2: As Fig. 1, but for the schematic IBM Hamiltonian of Eq. (9). 



In the physical region < x < 7r /2, H± exhibits a first order phase transition at 
Xc = arctanl/9 9 . As before, here we consider the interval — ir/2 < \ < 37r/2. 
The results for the distribution of ground state angular momenta are presented 
in Fig. 2. For N — 3k the ground state has L = in 75% of the cases and 
L = 2N in the remaining 25%. For N = 3k + 1 and N = 3k + 2 the ground 
state angular momentum is either L = (50%), L = 2 (25%) or L = 2N (25%). 
The variation in the L — and L = 2 percentages is due to the contribution 
of the (i-boson condensate, whereas the sum of the two is constant (75%). 

2.3 The SU{3)-SO{<6) case 

The transitional region between rotational and 7 unstable nuclei described by 
the Hamiltonian 



sin y 

+ 



(2s t -s t -d t -d t )(2S-S-d-d) 



N - 1 

+ (2s t xrf t ± v / 7d t xd t ) (2) -(2JxS±v / 7dxJ) (2) j , (10) 

does not show a phase transition in the physical region < x < tt/2 9 . Fig. 3 
shows that the distribution of the ground state angular momenta is very similar 
to the previous case. 
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Figure 3: As Fig. 1, but for the schematic IBM Hamiltonian of Eq. (10). 



3 Random interactions 

Finally, we apply the mean-field analysis to the general one- and two-body 
IBM Hamiltonian 



(11) 



in which the nine parameters of this Hamiltonian are taken as independent 
random numbers on a Gaussian distribution with zero mean and width a. The 
distribution of geometric shapes for this ensemble of Hamiltonians is deter- 
mined by the distribution of equilibrium configurations of the corresponding 
potential energy surfaces 

E(a) = <i4 sin a + 03 sin 3 a cos a + 02 sin 2 a + ao . (12) 

The coefficients a,i are linear combinations of the Hamiltonian parameters. The 
spectral properties of each Hamiltonian of the ensemble of random one- and 
two-body interactions are analyzed by exact numerical diagonalization 2 and 
by mean-field analysis 7 . 

In Fig. 4 we compare the percentages of L — and L = 2 ground states 
obtained exactly (solid lines) and in the mean- field analysis (dashed lines). 
There is a dominance of ground states with L = for ~ 63 — 77% of the cases. 
For N = 3k we see an enhancement for L = and a corresponding decrease 
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Figure 4: As Fig. 1, but for the random IBM Hamiltonian of Eq. (11). 

for L = 2. Also in this case, the equilibrium configurations can be divided into 
three different classes: an s-boson or spherical condensate, a deformed con- 
densate, and a (i-boson condensate. For the spherical and deformed solutions 
the ground state has L = (~ 63%) or L = 2N (~ 13%). The analysis of 
the ci-boson condensate is a bit more complicated due to the presence of two 
moments of inertia, Z5 and Z3. There is a constant contribution to L = 2N 
ground states (~ 10%), whereas the L = and L = 2 percentages show oscil- 
lations with N 7 . Just as for the schematic Hamiltonians, the sum of the L = 
and L = 2 percentages is constant and independent of N. 

4 Summary and conclusions 

In this contribution, we have investigated the origin of the regular features that 
have been observed in numerical studies of the IBM with random interactions, 
in particular the dominance of ground states with L = 0. 

In a mean-field analysis, it was found that different regions of the param- 
eter space can be associated with particular intrinsic vibrational states, which 
in turn correspond to definite geometric shapes: a spherical shape, a deformed 
shape or a condensate of quadrupole bosons. An analysis of the angular mo- 
mentum content of each one of the corresponding condensates combined with 
the sign of the relevant moments of inertia, provides an explanation for the 
distribution of ground state angular momenta of both schematic and random 
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forms of the IBM Hamiltonian. 

In summary, the present results show that mean-field methods provide 
a clear and transparent interpretation of the regular features that have been 
obtained before in numerical studies of the IBM with random interactions. The 
same conclusions hold for the vibron model 11 . For the nuclear shell model the 
situation is less clear. Despite the large number of studies that have been 
carried out to explain and/or further explore the properties of random nuclei 
no definite answer is yet available 12 . 
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